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Abstract: We compute the one- loop amplitude for a Higgs boson, a quark-antiquark pair 
and a pair of gluons of negative helicity, i.e. for the next-to-maximally helicity violating 
(NMHV) case, ^(ff, 1^, 2+, 3~, 4~). The calculation is performed using an effective La- 
grangian which is valid in the limit of very large top quark mass. As a result of this paper 
all amplitudes for the transition of a Higgs boson into 4 partons are now known analytically 
at one-loop order. 
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1. Introduction 

The hunt for the standard model Higgs boson is about to enter the endgame phase. The 
lower and upper limits coming from direct searches at LEP [1] and indirect constraints 
from precision electroweak data from the Tevatron and LEP [2] are now supplemented by 
the first direct limit from a hadron collider [3, 4]. With the increase of luminosity at the 
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Tevatron and the advent of running at the LHC, a discovery or a more stringent set of 
hmits is to be expected. 

An important search channel for the Higgs boson, in the mass range 115 < mn < 
160 GeV, is production via weak boson fusion [5]. A Higgs boson produced in this channel 
is expected to be produced relatively centrally, in association with two hard forward jets. 
These striking kinematic features are expected to enable a search for such events despite 
the otherwise overwhelming QCD backgrounds. Confidence in the theoretical prediction 
for the Higgs signal process is based upon knowledge of next-to-leading order corrections 
in both QCD [6-8] and in the electroweak sector [9, 10]. 

However, in addition to the weak process, a significant number of such events may 
also be produced via the strong interaction. In order to accurately predict the signal and, 
in particular, to simulate faithfully the expected significance in a given Higgs model, a 
fully differential NLO calculation of QCD production of a Higgs and two hard jets is also 
required. 

In the Standard Model the Higgs couples to two gluons via a top-quark loop. Calcu- 
lations which involve the full dependence on mt are difficult and a drastic simplification 

can be achieved if one works in an effective theory in which the mass of the top quark 
is large [11-13]. For inclusive Higgs production this approximation is valid over a wide 
range of Higgs masses and, for processes with additional jets, the approximation is justi- 
fied provided that the transverse momentum of each jet is smaller than nit [14]. Tree- level 
calculations have been performed in both the large-m^ limit [15, 16] and with the exact-m^ 
[14] dependence. 

Results for the one-loop corrections to all of the Higgs -|- 4 parton processes have 
been published in 2005 [17]. Although analytic results were provided for the Higgs qqqq 
processes, the bulk of this calculation was performed using a semi-numerical method. In 
this approach the loop integrals were calculated analytically whereas the coefficients with 
which they appear in the loop amplitudes were computed numerically using a recursive 
method. Although some phenomenology was performed using this calculation [18], the im- 
plementation of fully analytic formulae will lead to a faster code and permit more extensive 
phenomenological investigations. 

In recent years enormous progress has been made in solving the problem of evaluating 
virtual corrections to NLO scattering processes. Building upon the remarkable work of 
Bern, Dixon, Dunbar and Kosower during the mid-nineties [19, 20], unitarity construc- 
tions for these virtual corrections have developed into an efficient algebraic technique. The 
modern generalised unitarity method utilises quadruple cuts with complex momenta to 
freeze four dimensional loop momenta and uniquely determine the box coefficients [22]. 
The computation of triangle and bubble coefficients is also reduced to an algebraic proce- 
dure by application of an OPP style integrand reduction [23-25] or using direct analytic 
extraction [26] . Further developments employing D-dimensional cutting techniques [27-30] 
extend the method to compute full one-loop amplitudes. The procedure is well suited to 
numerical implementations and a number of automated approaches have been developed 
to the point of phenomenological applications [31-39]. 

In this paper we derive a compact analytic formula for the Higgs NMHV amplitude 
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with a quark-antiquark pair and two like-helicity gluons. This is achieved by sphtting the 
real Higgs scalar into two complex scalars {(f) and <p^) such that the Higgs amplitude is 
recovered in the sum [40]. Coefficients of box, triangle and bubble integrals are computed 
by applying the generalised unitarity method in four-dimensions. The rational terms are 
extracted from a Feynman diagram computation which is simplified using the knowledge 
of unphysical singularities in the cut-constructible terms. 

The paper is organised as follows. In Section ^ we describe the large top-mass ap- 
proximation and the decomposition of the Higgs into self-dual {(/)) and anti-self-dual {cj)^) 
components. Section |3| recalls the colour decomposition into primitive amplitudes and sec- 
tion ^ provides a guide to the current literature on Higgs + 4 parton amplitudes and recalls 
the known analytic results for 4'qqgg amplitudes that are needed to construct results for 
the Higgs boson. In section |5| we present analytic results that are sufficient for a complete 
description of the NMHV amplitude. We numerically evaluate the obtained expressions 
for the one-loop colour ordered amplitudes in section ^ before drawing our conclusions. 



2. Effective Lagrangian 



Our calculation is performed using an effective Lagrangian to express the coupling of gluons 
to the Higgs field [11], 

r 

(2.1) 



c 

— tr . 

This Lagrangian is obtained by replacing the full one-loop coupling of the Higgs boson to 
the gluons via an intermediate top quark loop, by an effective local operator. The effective 
Lagrangian approximation is valid in the limit mn < 2mt- At the order required in this 
paper, the coefficient C is given by [12, 13], 



C 



as 
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(2,2) 



Here v is the vacuum expectation value of the Higgs field, v = 246 GeV. The trace in 
Eq. (2T) is over the colour degrees of freedom which, since SU(3) generators in the funda- 
mental representation are normalised such that ti T'^T^ = 5°'^, implies that trG^^G'^'^ = 



Following reference [40] we will introduce a complex scalar field. 



\{H + iA), 



(2.3) 



so that the effective Lagrangian, Eq. (pT|), can be written as 

C 



£int 
H,A 



2 

c 



H tr Gf,^ Gf"" + iA tr G^^ *G'"' 



trGsn,.GZ + <l>^ trG 



ASD liv ^ ASD 



4 q T 



(2.4) 



where the gluon field strength has been separated into a self-dual and an anti-self-dual 
component. 



1 



{G 



fJ,U 



G 



pa 



(2.5) 



-3- 



Calculations performed in terms of the field (j) are simpler than the calculations for the 
Higgs boson and, moreover, the amplitudes for (/>^ can be obtained by parity. In the final 
stage, the full Higgs boson amplitudes are then written as a combination of (j) and 4>^ 
components: 



A{H, {pu}) = A{4>, {pk}) + A{4>\ {pk}) , 
A{A, {pk}) = -I {p,,}) - A{<^\ {pkS) 



(2.6) 



3. Definition of colour ordered ampfitudes 

The colour decomposition of the Hqqgg amplitudes is exactly the same as for the case 
qqgg which was written down in ref. [21]. For the tree graph there are two colour stripped 
amplitudes, 



Af{<t>,^,2,,2,,Aa) = Cg^ (T'^-(^)r'^^W)J^ 4 ^(</>, I5, 2„ C7(3), (t(4)) 
At one-loop level the colour decomposition is, 



(3.1) 



Nc J2 (T'^-(^'T«<^'*0j'^4;i(</',l,-,2„(T(3),a(4)) 

fTg52 



+ <5»3a4^Ji^^.3(^^l_^2,;33,43) 



(3.2) 



In these equations g is the strong coupling constant and cr is the ubiquitous one-loop 
factor, 

_ 1 r(i + 6)r^(i-e) 

= (4vr)2- r(l - 26) • ^^-^^ 

The colour stripped amplitudes ^44-1 and A^-^ can further be decomposed into primitive 
amplitudes. 



1 



A4;l{4>Aq,'^q,^gAg) = A^{(l),lg,2g,3g,4g) - ^Af{^,lq,2q,3g,4g 



+ ^Al{(P,lg,2„3g,ig), 



(3.4) 



and, 



A4;3('^, Ig, 3g,4g) = ^(0, I5, 2„ 3gAg) + Af{(f>, I5, 2g, 3g,4g) + A^ {cf> , Ig, 3g, 2g,4g) 

+ 4((/>, Ig, 2g, ig,3g) + ((/>, Ig, 2g,4g, 3 g) + A^ (cj) , Ig, 4g, 2g, 3 g) . 

(3.5) 



All of these colour decomposition equations, namely Eqs. (3A, |3.2| , 3.4 , |3.5|) are equally 
valid if the (f> is replaced by a (f>^ or a Higgs boson H. Sample diagrams contributing to 
each of the primitive amplitudes are shown in Figure ||. 
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(a) (b) 
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1 




(c) 






(d) 



Figure 1: Sample diagrams contributing to the primitive amplitudes (a) (0, 1^, 2^, 3g, 4g), (b) 
Af{(l),lq,2g,3g,4g), (c) {(/) , 1 g , 2 g , 3 q , 4 g) , Siud (d) ( (/) , 1 , , 2 , , 3g , 4g ) . The field can attach to 
any gluon line in the diagram. 

4. Known analytic results for Higgs + 4 parton amplitudes 

In this section we review results from the literature and collect formulae, for both tree and 
one-loop results, that will be useful in constructing the Higgs NMHV amplitude. 



4.1 Tree graph results 



The results for the tree graphs that are primarily of interest here, i.e. (pqqgg amplitudes 
with gluons of the same helicity, are: 



-lA 



(°)/'^ ir 2+ 3~ 4" 



(3|(1 + 4)|2]^(41) 

[2 4] S124 

(4|(l + 3)|2f (13) 
[2 3] S12S123 



1 1 

S12 S41 



+ 



(l|(3 + 4)i2]^ 
(12) [2 4] [2 3] [34] 



-iA 



ir 2+ 3+ 4+ 

'^j ^q 1 '^q J '^g 1 ^g 



and for the subleading colour piece. 



-iA 



(0), 



-iA 



1 r 9^ .S+ 4~ 



^> 1 2g , 3g , 4^ 



(4|(l + 2)[3]^ (2|(1 + 4)|3]^ 
[12][2 3]si23 [34][41]s34i 



0. 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



A brief summary of our spinor notation is given in Appendix Compact analytic expres- 
sions for all helicity amplitudes are presented in references [45, 46]. 
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By using parity and charge conjugation [47], we can relate these (jjqqgg amphtudes to 
ones for (p^qqgg with the same hehcity assignments of quark and antiquark. This relation, 
valid at any order of perturbation theory, n, reads, 



4")(,/.,2^'^^4^4;'^^3;'^3^ 



(4.5) 



We thus see that the cj)^ amplitude in which we are interested is zero. 



-iA 



(0), 



5 ) 2g , 3g , 4g ) — 



(4.6) 



so that, at tree graph level, the NMHV Higgs amplitude in which we will ultimately be 
interested is simply given by Eq. ( [4.11 ). 



4.2 One-loop results 

We begin this section with a brief survey of the literature, to indicate where original 
calculations of Hqqqq and Hgggg amplitudes may be found. We then turn to the Hqqgg 
amplitudes, concluding the section by quoting the results for the (f^ amplitude that must 
be combined with the new (j) amplitude calculation that we present here. 



4.2.1 Hqqqq amplitudes 

The full one-loop results for this process, both for pairs of identical and non-identical 
quarks, are already available in the literature. The matrix element squared has been 
computed in ref. [17], with results for the amplitude presented in ref. [46]. 



4.2.2 Hgggg amplitudes 

In principle there are 16 combinations of amplitudes, but this number is reduced to four 
independent amplitudes by parity and cyclicity. The references to the complete set of 
needed amplitudes are given in Table ||. In addition a nice summary of all the one-loop 



H amplitude 


4> amplitude 


(f)^ amplitude 


^(i/,l+,2+,3+,4+) 
^(//,l,-,2+,3+,4+) 
^(i?,l-,2-,3+,4+) 
^(i7,l-,2+,3-,4+) 


^(0,1+2+, 3+4+) [44] 
^(0,1^,2+, 3+ 4+) [44] 
^(0,1,, 2-, 3+4+) [42] 
^(0,l-,2+,3-,4+) [43] 


^(0t,i+ 2+,3+ 4+) [41] 
^(</<t,i-,2+,3+ 4+) [45] 
^(0t,l-,2-,3+ 4+) [42] 
^(</<t,l-,2+,3-,4+) [43] 



Table 1: (f) and (f)^ amplitudes needed to construct a given one-loop Hgggg amplitude, together 
with the references where they can be obtained. In all cases the 0^ amplitudes are constructed 
from the amplitudes given in the reference using the parity operation. Results for all helicity 
combinations are also written, in uniform notation, in ref. [45]. 



results for the Higgs -|- 4 gluon amplitudes is given in ref. [45]. 
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H amplitude 


(j) amplitude 


cj)^ amplitude 


^(i/,ir,2+ 3+,4+) 
^(i/,ir,2+ 3-,4-) 
A{H, ir,2+ 3+ 4-) 
^(i7,l,-,2+ 3-,4+) 


^(0,1^,2+, 3+ 4+) [44, 46] 
^(0, ir,2+,3-,4-) 
^(0,15,2+, 3+ 4-) [46] 
^((^,15,2+, 3-, 4+) [46] 


^((/>t,l^r,2+,3+,4+) 
^((/>t, 15,2+, 3-, 4-) [44, 46] 
^(<At,l,-,2+,3+,4-) [46] 
^(</.t,ir,2+,3-,4+) [46] 



Table 2: <j) and 0^ amplitudes needed to construct a given one-loop Hqqgg amplitude, together 
with the references where they can be obtained. In all cases the (f^ amplitudes are constructed from 
the (p amplitudes given in the reference, using the parity operation. The cases where the gluons 
have the same helicity, which have no associated references, are the subject of this paper. 



4.2.3 Hqqgg amplitudes 

In principle there are 8 combinations of amplitudes, since helicity is conserved on the quark 
line, but because of parity invariance only four Higgs amplitudes are independent. The 
references to the amplitudes already calculated in the literature are given in Table |. From 
this table we see that the Higgs amplitude A{H, 1^, 2+, 3~, 4~) requires, in addition to the 
calculation of a previously unknown (j) amplitude, also the results for the corresponding (j^ 
amplitude from ref. [46]. 

The (j)^ results that we shall need can be derived from the following amplitudes in the 
case of A4.1, 



ir,2+,3+,4+) = 2iAf (</.t,ir,2+,3+,4+) 



1 

+ 2 



1|(2 + 3)|4] (12) [2 3] (31) 



(2 3) (3 4) 
1- 2+ 3+ 4+'> = 



1~ 2+ 3"*" 4"'"'> 



+ 
1 

2 



(2 3) (3 4) (41) 



(1|(2 + 3)|4] 
(2 3) (3 4) 
1(13) [3 4] (41) 
3 



+ 



1 (13) [3 4] (41) 
3 (12) (3 4)2 
(12) [2 3] (31) 
(2 3) (3 4) (41) 



(12) (3 4)2 ' 

whilst the subleading partial amplitude A4.3 also requires the results. 



(4.7) 
(4.8) 
(4.9) 



-af(,/.,ir,2+ 3+,4+) 
4(,/.,ir,2+,3+,4+) 



(1|(2 + 3)|4]2 1 
(12)(2 3)si23 2 



(13)(1[(3 + 4)|2] 
(2 3) (3 4) (41) 



+ 



0. 



(13)^ [34] 
;i2) (2 3) (34)_ 

(4.10) 



To obtain the form that is most useful for the calculation of A{H, 1^, 2+, 3~, 4^"), we 
relate the (p^qqgg amplitudes to the 4>qqgg ones by using the relation in Eq. (4^). Thus 
we obtain the required results by performing the transformation 1 <-> 2, 3 ^ 4, () ^ [] and 
reversing the sign. The amplitudes contributing to A^-i are, 



1 ~ 2+ 3~ 4~ 

^^q T^g ^ ^g 



2iA 



(0), 



1-0+0 



4- 
9 ' 9 



1 

+ 2 



(3|(1+4)|2] , [21] (14) [24] 



+ 



[14] [3 4] [14] [3 4] [2 3] 



1 [2 4] (34) [2 3] 
3 [12][34]2 



-af(</>t,ir,2+,3,-,4-) 



(3|(1+4)|2] , [21] (14) [24] 



+ 



[14] [3 4] [14] [3 4] [2 3] 



(4.11) 
(4.12) 
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iA{{4>\l^X,?,-A-„ 



1 [2 4] (34) [2 3] 



^<?'-</'"9-9^ 3 [12][34]2 ' 
while the additional subleading contributions become, 



(4.13) 



Mi(</.t,ir,3-,2+ 4-) = 2 



(3|(4 + 1)|2]^ 1 
[2 4][41]si24 2 



[21](4|(1 + 3)|2] , [12]2(13) 



+ 



[41] [13] [3 2] [2 4] [41] [13] 



= i A^^^ {(j), 1^ ,3g ,2^ ,4g) + terms antisymmetric in {3^4}. 

(4.14) 

We note that all of these amplitudes are finite because of the vanishing of the corresponding 
tree graph results (see section . 



5. One-loop results 

In this section we present analytic expressions for the full one-loop corrections to the 
process ^4^^((/), 1^, 2+, 3~, 4~). All expressions are presented un-renormalised in the four- 
dimensional helicity (FDH) scheme (setting = 0) or 't Hooft-Veltman scheme (setting 
Sr = 1). 

We employ the generalised unitarity method [22, 48-51] to calculate the cut-constructible 
parts of the left-moving, right-moving and nj one-loop amplitudes. This relies on the ex- 
pansion of a one-loop amplitude in terms of scalar basis integrals, 

^cut-cons. 1_ ^ 2+ 3^-, 4p = ^ C,,a4■,^ + ^ ^3;.:^3;^ + Yl ^^^^^^^ (^'l) 

i i i 

In this sum each j-point scalar basis integral {^j-i) appears with a coefficient Cj-^i. The 
sum over i represents the sum over the partitions of the external momenta over the j 
legs of the basis integral. Multiple cuts isolate different integral functions and allow the 
construction of a linear system of equations from which the coefficients can be extracted. 
We use the quadruple cut method [22] which freezes the loop momenta and determines each 
box coefficient uniquely. Triangle coefficients are determined using the Laurent expansion 
method [51], whilst the two-point coefficients are determined via Stokes' Theorem applied 
to functions of two complex-conjugated variables [52]. Results were obtained using the 
QGRAF [53], FORM [54] and S@M [55] packages in order to control the extensive algebra. 

5.1 Results for ^4;i((/), Ig, 2g, 3^", 4^ ) 

The partial amplitude A4^-i{(j)^ lg,2g,3~,4~) is calculated from three primitive amplitudes 



according to Eq. (3^). We shall deal with each of these ingredients in turn. 



5.1.1 ^f(,^, 1^,2+, 3-, 4-) 

The full result for this primitive amplitude is given by, 
- iAl (</), ir , 2+ , 3^- , 4- ) = -iAf ir , 2+ , 3" , 4; ) x V^- 
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+ 



'134 



[14] [34] (2|(1 + 4)|3] 
(1|(3 + 4)|2]2 



2mh . 



Ls_i(si4, S34; S134) + Ls_i (si2, S134; S34, m; 



(1|(2 + 3)|4][23] [34] 



2mh , 



Ls_i(s34, S23; S234) + Ls_i (si2, S234; S34, rrij^) 



+ 



+ 



+ 



ml(14)2(2 4) 



(3|(1 + 4)|2]-^ 



(12) (2|(1 + 4)|3](4|(1 + 2)|3]si24 [12] [24] (3|(1 + 2)|4]si24 



[23]^(4|(2 + 3)|1]^ 



[12][13]3(4|(l + 2)|3]si23 (12)(l|(2 + 3)|4](3|(l + 2)|4]si23 



Ls_i(si2,Sl4;Sl24) 
Ls_i(si2,S23;Sl23) 



(4|(1 + 3)|2]^ 



[12][2 3](4|(l + 2)|3]si23 (12)(l|(2 + 3)l4](3|(l + 2)|4]si23 



'2mh , 



'2mh , 



Ls_i (S34,si23;si2,"2j,) +Ls_i (si4, ^123; S23, "I. 



+ 



^^(14)' (2 4) (3|(1 + 4)|2]'(3|(2 + 4)|1] 

(12) (2|(1 + 4)|3](4|(1 + 2)|3]si24 [12] [14] (3|(1 + 2)|4]si24 



'2m/i , 



'2mh , 



Ls_i (s34,Sl24;Sl2,rn^) +Ls_i {S23,si24;si4,mj 



a 



3;0|12|34l 



ir , 2+ , 3- , 4- (si2 , S34 , ) - ^3.^141123 (0, 1^ , 2+ 3; , 4; jl^- (^23 ,^14, 



2(13)2(34) (4|(1 + 2)|3] [12] 



L3(S123,S12) 



, (34) (3 1) ((4|(1 + 3)|2] [1 3] - 3(4|(2 + 3)|1] [23]) ^ 
T TTTTrrr. '-'2{Si23, S12) 



6 [31] 



+ 



(1 3) (l6(4|(l + 3)|2]'[1 3]' - 3(4|(1 + 3)|2](4|(1 + 2)|3] [2 1] [3 1] + 6(4|(2 + 3)|1]'[2 3]') 



6si23[3 1]'[3 2] 



XLi(si23,Sl2) 

2si24{34)'(14) [42] 



L3(S124,S12) 



- (34) (14) »±iMi4^^ffll±Wl£,(,,,,,,.,, 



+ 



6 [41] 

(3|(1 + 4)|2](9 S124 (34) [2 1] - 22 (3|(1 +4)|2] (4 2) [12]) 

6si24[41] [21] 
(14) (13) (4|(2 + 3)|1] [12], 



Ll(si24, S12) 



2 [31] 



-L2(S123,S23) 



(13) (41(2 + 3)11, 3 (4|(1+3)|21|13|+2(4|(2 + 3)|11[23| ^_^^^^^^^^^^ 

2S123[13] 



2 [4 3] 

+i?^(0,ir,2+ 3;,4;), 



2 [4 3] 



(5.2) 



and the remaining rational terms given by. 

(34)(3|1 + 4|2](2(24)[42]-(12)[21] 

i?^(</.,l^,2+ 3-,4 ■ 



^' ^' 12^124 (12) [21] [41] 

(2 3) (4|(1 + 3)|2]2(3(12) [21] - 2(2 3) [3 2]) - 2(13)^(2 4) (4[(2 + 3)|l] [21] [3 2] 



12.si23 (12) (2 3) [2 1] [3 1] [3 2] 
, 5(3 4)^ ^ 5(34) (4|(1 + 3)|2] ^ (4|(1 + 3)j2]2 



12 (2 3) [31] 6(2 3) [31] [3 2] 6 (1 2) [2 1] [3 1] [3 2] 

(13) (14) (2 4) [21] (13) (3 4) (3 4)^ [4 2] ^ (13) (24) (4|(1 + 3)|4] 



3(12) (2 3) [31] [3 2] 12 (12) [41] 6 (1 2) [2 1] [4 1] 4 (1 2) (2 3) [3 1] [4 3] 
(13)(4|(1 + 3)|4] 5(14)^[41] (14)^[42] 

3(12) [41] [43] 12(12) [31] [43] 6(12) [3 2] [43] ■ ^ ' 

The coefficients of the three mass triangles were calculated using the method of ref. [26] , 

_ + _ _ _ ml{2,A)\lK\f 
C3;,|i,34(0, 1, , 2, , 3, , 4, ) - ^(^_^p^i2)(3E:J>(4Kt> ' ^'-'^ 

with Ki = —pi — P2 — P3 — Pi, K2 = —pi — P2 and the massless vector K\ given by, 
and where 7 is given by the two solutions, 



7± = i^i • ± V • ^2)2 - KlKl . (5.7) 
The other triangle coefficient is, 

_ + _ _ _ ml{lAf{'iK\f 
C3;,|4i|23(0, 1, , 2, , 3, , 4 J - - 27(7-m^)(li^J>(2i^J) ' ^'''^ 

with Ki = —pi — P2 — — Pi, K2 = —pi — Pi and K\ given in terms of these vectors by 
Eq. (P). 



~2m/i 

The definitions of the box integral functions Ls_i and Ls_]^ can be found in Ap- 
pendix 1^, together with expressions for Li, L2 and L3. In addition to logarithms and 
polynomial denominators, the latter functions also contain rational terms that protect 
them from unphysical singularities. Thus, for example, 

^ , , log(s/t) 1 (\ 1\ 

which is finite in the limit that s ^ t. 
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5.1.2 A^{^,l^,2+,3g,4g) 

The result for the right-moving amphtude, Af{(j), 1^,2^,3~,4~) is, 

-af(0,l^r,2+ 3^-,4;) = -af (c/>,l^r,2+ 3^-,4^-) x y« 
, [12]^(4|(1 + 2)|3]\ , , (3|(1 + 4)|2]\ 

+ [13]^[23].,,3 ^ [14] [2 4] .,.4 ^"^^ 

(1|(3 + 4)|2]2 ~2m/^ 2x 

Ls„i (Sl4,S234;S23,'W<Aj 



[23][34](1|(2 + 3)|4] 

2 

^134 ~2m/i 2x 

+ [14][34](2|(1 + 4)|3] ^'-1 (^23,^i34;.i4,m<^) 

_ , _ _ 3^ 2n 1 (14)^[12]^(3|(1 + 2)|4]2 ^ 

- C'3;(/.|41|23('A, Ig , 2g , 3g , 4g jig (S23, ^14, - ^ [1 4] [2 4] gi24 ^2(5124, SI2) 

g (34)(3|(l + 4)|2] ^ 1 (31(1 + 4)12]^ ^ 

+ ^ P '^^^ + 2 [14] [2 4] .124 °^ '''' 

1 (14)^[24]^.i3, ^ (34) (1|(3 + 4)|2] ^ 

" 2 [23] [34] (1|(2 + 3)|4] ^2(^234, .23) - 2 Li(s234,.23) 

1 (1|(3 + 4)|2]2 ^ i((12)[12](4|(2 + 3)|l])'[23] ^ 
+ 2 [23] [34] (1|(2 + 3)|4] ^0(^234,. 23) - ^ ^^j3^^^^ L2(.i23,523) 



+ 



(13) [12] (4|(1 + 2)|3](4| (2 + 3)[1] 
(2 3) [13]^ [2 3] 

^ (13)[12](4|(1 + 2)|3](4|(2 + 3)|1] ^ 1 (4|(2 + 3)|1]^ [23] 

si23[13]'(2 3) [2 3] 2 [13]3si23 

1^(13) [12] (4|(1+2)|3])\ 

2 [13] [2 3]S123 L2(S123,S12) 

-2( 13) [13] - (2 3) [2 3]) 
(2 3) [13]^ [2 3] 



Lo(si23,S23) 



+ (34) [12] (4[(1 + 2)|3] ^--^-:1:.72L:. U{s^2S,s,,) 



^ h 91 Min ^ oM.i (23) (4|(1 + 3)|2]+2(13) (4|(2 + 3)|1] ^ 

+ [12] (4|(1 + 2)|3] Lo(si23,si2) 

[13] (2 3)[2 3]si23 

+ i?^(0,l^,2+,3^-,4g) , (5.10) 
with 

e2 V-si2y 2e V-si2y 2 2 ^ ^ 

The remaining rational pieces in Eq. ( ^.lOD have the following form: 

R-i^, l,-,2-,3-,4,-) = _MW ^ (4|(1 + 2)|3F[21]- _ (14)^ [21] 
^ " '^^^^ 2(23)[31]^ 2si23[31]='[3 2] 2 (1 2) [3 1] [3 2] 

[21] ((13)2 (23) (4|(1 + 2)|3]2 [31]' + (12)=^ (4|(2 + 3)|1]2 [21] [3 2] 
^ 45f23 (12) (2 3) [31]=^ [3 2] 
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+ 



2si24[41][4 2] 2(12) [41] [4 2] 
(13) (14) [4 2] 5234 (1 4)2 [4 2]2 



+ 



(3|(1 + 4)|2]2 (13)2 [21] , ( 14)2(31(1 + 2)14]^ [21] 

4.22^(12) [41] [4 2] 
(14)2 [4 2]2 



2(1|(2 + 3)|4][4 3] 4(2 3)(1|(2 + 3)|4][3 2]2[43] 2(1|(2 + 3)|4] [3 2] [4 3] 



(5.12) 



5.1.3 A 



1 ~ 2+ 3~ 4" 



MvV) ^q : "q J '-'g : 

The fermion loop contribution is, 



-iA\ 



^ Ig ) 2g , 3^ , 4g 



\(t>, Ig ,2^ ,3g , 4g 



2 
3e 



-Sl2 



10 



+ ^(13)2 (34) [12] (4|(1 + 2)|3] £3(5123,^12) + ^(14)2 (34) [12] (3|(1 + 2)14] £3(5124,512) 

+ ^ (13) (34) (4|(1 + 3)|2] £2(5123,512) + ^ (14) (34) (3|(1 + 4)|2] £2(5124,512) 

2 (13) (4|(1 + 3)|2]2 ^ ^ , , 2 (14)(3|(1 + 4)|2]2 ^ ^ 

" 3 (12) [2 3] [12] Li (5123, 512) + 3 (12) [2 4] [12] ^1(5124, 512) 



+ 



2 (13) (4|(1 + 3)|2]2 

3 (12) [12] [2 3]si23 
(13) (34) (4|(1 + 3)|2] 



Lo(5l23,5l2) + 



6(12) [12]s 



123 



2 (5i2 + 5i4)(3|(l+4)|2]- 

3 (12)[14][2 4][12]si24 
14) (34) (3|(1 + 4)|2] _ 1 (13) (14) 
6(12) [12]si24 3 (12) [34] 



Lo(5l24,5l2) 



(5.13) 



5.1.4 Relation for rational terms 

We note that the rational terms in the three leading colour primitive amplitudes obey, 

n[A^{<f>, Ig, 2g, 3g,Ag) + Af{(f,, Ig, 2g, 3g,4g) + 4('/'> 3g,4g)^ 

+ 2Af\^\lg,2g,3gAg)=0, (S.M) 

a formula analogous to that found in super-symmetric decompositions of QCD ampli- 
tudes [21]. This property is helicity independent and has also been checked for the previ- 
ously known MHV amplitudes [46]. For the NMHV helicity assignment at hand, namely 
(1^, 2+, 3g , 4~), we note that the tree graph result that appears in Eq. ( ^.14|) is zero (c.f. 
Eq. ([4.6|)). We stress that the TZ operation extracts the full rational term, including com- 
pletion terms from the functions L3 and L2. Thus it corresponds to dropping all logarithms, 
box functions and TZ-functions. 

We conclude this section by noting that the three primitive amplitudes for the helicity 
assignment (1^, 2+, 3+, 4+) displayed in Eqs. (|3), (|]|) and (|]|), also satisfy Eq. (|5.14| ). 



For these amplitudes, which are purely rational, the TZ operation leaves the amplitude 
unchanged. 



5.2 Results for A^-i 



S Ig, 2g, 3g , Ag 



We can calculate the result for ^44^3 using Eq. ( |3.5D . Given the results for Aj^ and Aj^ in 
the previous section the only missing ingredient is A^{cl), 1^, 2~, 3^, 4~). 
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5.2.1 Box-related terms for A^((j),lg ,2~ ,3^ 

The calculation of the box-related terms in (f)qgqg ( 1 — ) is easily performed using the 

methods given in ref. [22]. The result is, 

-a^''^°"(<^,ir,2;,3+ 4^) = -iAf\ct>,l^,2;,3+,4;)xVi^ 
(12)H 



+ 
+ 

+ 



(l|(2 + 3)|4](3|(l+2)|4]si23 
[34]2(1|(3 + 4)|212 



(2|(1+4)|3]\ 

Ls_i(Si2,S23;Sl23) - r-, .1 ro .1 Ls_i (Si4, S34; S134) 



[14] [34]S134' 



[32] [24]^(1|(2 + 3)|4] 
(3|(2 + 4)|1]. 



Ls_i(s23,S34;S234) + 



s2 



[12] [24](3|(1 + 2)|4] 



Ls_i(si2,Si4;Si24) 



2 

124 



'2mh 



[12] [14](3|(1 + 2)|4](3|(1+4)|2] 



Ls_i (s23,Sl24;Sl4,mj 



(11(2 + 4)13] 



'2m/i 



+ 



[23] [34](1|(2 + 3)|4](1|(3 + 4)|2] 
ml{12f 



S123 



+ 



Ls_i (si2,S234;S34,mj,) 



(4|(1 + 2)|3]^ 



(3|(1 + 2)|4](1|(2 + 3)|4] [12][23] 



— 2mh 2 \ ' — 2m/i q s 

Ls_i (S34,Sl23;Sl2,m^) +Ls_i {su,Si23;S23,m^) 



(5.15) 



with 



2 \ e 



-S34 



2 \ en 



-S41 



1 ^ 



3e V-S123 



_l £i 

4 3 



(5.16) 



As we shall see in the next section, no further information is required for the calculation 
of the A4.3 which is completely determined by box diagrams alone. 



5.2.2 Full result for ^4;3 

The full result for the partial amplitude ^4;3 is, 

- «^4;3(0, l^,2+,3~,4~) = -zAf^(0, l^,2+,3~,4~) X ^5(512, S34, Sl3, S24) 



+ 



S123 



(4|(1 + 2)|3]2[12]2^ [2 3]2(4|(2 + 3)|1]3 



[13f [2 3] [13]^ [12] (4|(1 + 2)|3] (12) (1|(2 + 3)|4](3|(1 + 2)|4] 



X Ls_l(si2,S23;Sl23) 



+ 



+ 



Sl24 
1 

S123 



mj (14)2 (2 4) 



(3|(1+4)|2]^(3|(2 + 4)|1] 
(1 2) (2|(1 + 4)|3](4|(1 + 2)|3] [1 4] (3|(1 + 2)|4] [1 2] 



Ls_i(si2,Si4;Si24) 



(41(1 + 3)|2]^ 

(1|(2 + 3)|4](2|(1+3)|4] [13] [23] 



Ls_i(si3,S23;Sl23) 



'341 



+ 



[13] [34] (2|(1 + 3)|4] 

^341 

[14] [34] (2|(1 + 4)|3] 



'2mh , 



Ls_i(si3,Si4;S34i) +Ls_i (S23,S341,S14,"7,<J 



2mh , 



Ls_i(si4,S34;S34i) + Ls_i (si2,S341,S34,m0) 
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(1|(3 + 4)|2]^ 



-2mh , 



+ 



[23][34](1|(2 + 3)|4] L 

[24f (1|(2 + 4)|3]^ 
[23] [34]3(1|(2 + 3)|4] 



Ls_i(s23, S34; S234) + Ls_i (si2, S234, S34, m^' 

Ls_i(s23,S24;S234) 



(l|(3 + 4)|2p 



-2mh 



(1|(2 + 4)|3] [24] [34] 



Ls_i (S14,S234,S23,"1. 



+ 



S123 



m;(13)^(2 3) 



(4|(1 + 3)|2]^ (4|(2 + 3)|1] 
(12)(2|(1 + 3)|4](3|(1+2)|4] ' (4|(1 + 2)|3] [12][13] 



+ 



-2111/1 



+ 



X Ls_i (S14,S123,S23,"1^) 

m^(13)3 



S123 



(4|(1 + 3)|2]= 



12)(3|(1 + 2)|4](1|(2 + 3)|4] (4|(1 + 2)|3] [1 2] [2 3] 

2m/i 



+ 



X Ls_i (s24,Sl23,Si3,m^ 

m4,(13)^ 



S123 



(2|(1 + 3)|4](1|(2 + 3)|4] 



+ 



(4|(1 + 3)|2]- 
[13] [2 3] 



Ls_i (S34,Sl23,Sl2,'n^w 



+ <^3^4 



(5.17) 



where the function containing poles and associated logarithms is conveniently written as, 



V5(S12, S34, Sl3, S24) 



1 



-S\2 



+ 



-S34 



-Sl3 



-S24 



(5.18) 



We note that the apparent double pole in e in Eq. (|5.18|) is cancelled upon expanding about 
e = 0. 

This result for the ^ amplitude is particularly simple, containing neither bubble contri- 
butions nor rational terms. This is also true for the helicity amplitude ^4.3(1;^, 1^, 2^, 3~ , 4+) ^, 
which can easily be checked using the previously calculated results in ref. [46]. It is there- 
fore more efficient to program the full result for ^4^3, rather than to program the individual 
primitive amplitudes using Eq. (|3.5|). 

Furthermore, for the case of two negative gluon helicities calculated here one can check 
using Eq. (|3.5|) and Eqs. ( 4.1l| ), (|4.12D , ( [4.14 ) that the corresponding <^'^ amplitude is zero. 
Therefore we have, 



^4;3 (/f, I5 , 2+ , 3- , 4- ) = i A4;3 (^,15,2+ 3" , 4" ) = A4;3 (</>, I5 , 2+ , 3" , 4^ 



(5.19) 



6. Numerical results 



Here we present evaluations of the new amplitudes at the same kinematic point as used 
previously in the literature [17, 46]: 



(-1.0000000000, 0.00000000000, 0.00000000000, 0.00000000000), 



""^The amplitude j44;3((^, l5,2q,3^,4g ) is not independent and is obtained by swapping labels 3 and 4. 
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(<A, 1^,2+, 3-, 4-) 


id^ 1- 2+ 3" 4~) 






1/e 






/l(O) 








+6.49907535901 

— Z.oyoUol44olD I 







-19.49722607702 
4-7 ^ 7Q9 AA'iAAA7 i 


-64.62496304875 
zl^ 7fi1 1 9071 ^71 -i 


-31.60558356648 
1 '^7 ^fifi'^Q'^ni 4^9 9 


— i / .ooo4yzuouuo 
+6.14361664194 i 


aR 

^4 


-6.49907535901 
-1-9 3930814481(17 


-23.12631834140 
— 1 4 67090390044 9 


-48.74190400225 
—39 06965559875? 


-1-4. ^8546771410 
-1.38718545292 i 


^4 





-4.33271690600 
+1.59538763210 i 


-14.98058321393 
-9.88874973495 i 


—0 9981 964091 9 
+0.02973170727 i 




-18.77510659269 
+6.91334640579 i 


-69.27656696526 
-41.47211867326 i 


-51.15745514499 
-149.70134751402 i 


-18.24519911293 
+6.34730120438 ? 


^4;3 





+2.61083477136 
-0.05119106396 i 


+17.75737443413 
+4.93097014463 i 






Table 3: Numerical values of (j)qqgg and (f>^qqgg primitive amplitudes (above) and the amplitudes 



multiplying the two different colour structures (below), at the kinematic point defined in Eq. (6.1) 



ki = (0.30674037867, -0.17738694693, -0.01664472021, -0.24969277974), 

A;2 = (0.34445032281, 0.14635282800,-0.10707762397, 0.29285022975), (6.1) 

h = (0.22091667641, 0.08911915938, 0.19733901856, 0.04380941793), 

ki = (0.12789262211, -0.05808504045, -0.07361667438, -0.08696686795). 

We have used a scale /i = mn, set 5r = 1 (corresponding to the 't Hooft-Veltman scheme) 
and, in assembling the amplitude A^-i, have used Uf = 5. The results for the final Higgs 
amplitudes presented in Table ^ agree with those from the semi-numerical calculation of 
ref. [17] to one part in 10®. Note that these results depend on an overall phase that can 
be removed by dividing out by the corresponding Born calculation. Using the analytic 
expressions for all the Hgggg, Hqqgg and Hqqq'q' amplitudes that are now available we 
can also confirm^ the numerical values for the matrix elements squared given in ref. [17]. 



7. Conclusions 

In this paper we have computed the last remaining, analytically unknown, helicity am- 
plitude contributing to the NLO corrections to Higgs plus two jet production at hadron 
colliders. This builds upon the previously known semi-numerical results [17] and completes 
the set of compact analytic formulae [41-46]. 

We employed a generalised unitarity approach to calculate the cut-constructible parts 
of the amplitude. Completion of the logarithmic terms to remove unphysical singularities 

^Fortran code that calculates all the amplitudes can be downloaded from mcfm.fnal.gov. 
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was used to simplify rational terms extracted from a Feynman diagram calculation. Sim- 
plifications in the construction of the subleading colour amplitude, A4.2, showed that this 
component is free from all rational, bubble and triangle terms. Similar relations between 
the rational terms in the left, right and fermion loop primitive amplitudes were used to 
find a compact analytic structure. 

Our results have been verified against the known numerical results and we envisage that 
they will provide the means for a faster and more flexible analysis of Higgs phenomenology 
at the Tevatron and the LHC. 
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A. Spinor notation 

Our spinor notation is quite standard in the QCD literature, (for a review see refs. [56, 57]) 

The function u±{ki) is a massless Weyl spinor of momentum fcj and positive or negative 
chirality. In terms of these solutions of the Dirac equation, the spinor products are defined 

by, 

{ij) = = = u-{ki)u+{kj) , (A.l) 

[ij] = = (^+lr ) = u+{ki)u^{kj) . (A.2) 

We use the convention [ij] = sgn{k^kj) {ji}*, so that, 

{ij) [j i] = ■ kj = Sij . (A.3) 

We further define, 

Siji = {ki + kj + kif = {i j) [j i] + {j I) [I j] + {i I) [I i] , (A.4) 

{a\i\b] = {ai)[ib], {a\{i + j)\b] = {ai) [ib] + {aj) \j b] , (A.5) 

{ji) = -{ij), [ji] = -[ij]- (A.6) 

Simplification of the formula can sometimes be achieved by using the Schouten identity. 



(a 6) {cd) = {ad) {cb) + {ac) {bd) , (A.7) 
[a b] [c d] = [a d] [c b] + [a c] [b d] . (A.8) 



B. Definitions of special functions 

With the definition of the Lj(s,t) basis functions, 



we can define the completions, 



L3(.,t)-L3(.,t) = L3(.,t)-^^^(i + i 

U{s,t) ^Li{s,t) = U{s,t) , (B.2) 

such that the completed functions arc finite in the limit that s ^ t. Note that the definition 
of these completed logarithmic functions is similar in spirit, but different in detail from the 
definitions in ref. [46]. 

We also need the box functions from the scalar box with one massive external leg, 

Ls_,(M;m^) = Li, (l - ^) + U, (l - ^) + log (^) log (;^) - ^ , (B.3) 

and coming from the box with two adjacent massive external legs, 

Ls'r,\s, t;mlml) = -Li, ( 1 - H^i) - LiJ 1 - ^) - I log' 



t V t 2 ° \-t 



1, / -s 



where the dilogarithm is defined as usual by, 

Jo y 

Ig"^ is the three mass triangle function defined, for example, in Eq. (II. 9) of ref. [58], 

1 

l3™(si2,S34,S56) = / d^tti 6(1 - tti - a, - as) 7. (B.6) 

Jo [-5120102 - 5340203 - 5560301] 

Explicit results for this integral can be found in refs. [59], [60] and [61]. 
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